Real-Time Imaging of K atoms on Graphite: Interactions and Diffusion 
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Scanning tunneling microscopy (STM) at liquid helium temperature is used to image potassium adsorbed on 
graphite at low coverage (^0.02 monolayer). Single atoms appear as protrusions on STM topographs. A sta- 
tistical analysis of the position of the atoms demonstrates repulsion between adsorbates, which is quantified by 
comparison with molecular dynamics simulations. This gives access to the dipole moment of a single adsorbate, 
found to be 10.5=bl Debye. Time lapse imaging shows that long range order is broken by thermally activated 
diffusion, with a 32 meV barrier to hopping between graphite lattice sites. 



The addition of alkali atoms onto and into graphite has long 
been recognized as a powerful technique for engineering its 
electronic properties. Beyond the simple doping that results 
from a charge transfer, alkali intercalation has been shown to 
lead to superconductivity in several graphite compounds d 
E). The emergence of graphene brought a renewed interest in 
the field: possibilities range from simply achieving very high 
doping levels Q, to so-called atomic collapse for adatoms that 
transfer sufficient charge to the substrate (4), to inducing two- 
dimensional superconductivity in doped graphene |5j|6]]. 

For all the recent excitement about alkali-on-graphite/ene 
systems, little is known about interactions between the 
adatoms and the surface, or between the adatoms themselves. 
These are questions that must be resolved to better understand 
the effects of alkali doping. Coulomb interactions between 
adatoms are influenced by screening in the graphite, which it- 
self is a topic of great interest as it relates to electron-electron 
interactions in this system. An understanding of adatom in- 
teractions is also important to the study of phase transitions in 
two dimensional systems, and of adsorbate self-organization 
phenomena on surfaces EE). 

The structural phase diagram of alkali adatoms on graphite 
has previously been explored using low energy electron 
diffraction (LEED)[9|. For the case of potassium, a coverage 
higher than 0.1 monolayer (ML) leads to the condensation of 
a metallic 2x2 phase. In the dilute limit more commonly used 
for electronic doping of surfaces, K adatoms form a dispersed 
phase, with a large spacing that is believed to result from elec- 
trostatic repulsion between the partially-ionized adatoms af- 
ter charge transfer to the graphite l9l. Although spectroscopic 
and structural studies (LEED, Electron Energy Loss Spec- 
troscopy, Photoemission) support this general picture I9HTT1I. 
local probes have the potential to offer much more detailed 
insight into adatom interactions [12]. 

In this letter we take advantage of the real space imaging 
capability of the scanning tunneling microscope (STM) to ex- 
plore the spatial distribution and dynamics of potassium ad- 
sorbed on graphite in the dilute limit. Images of single station- 
ary K atoms on highly oriented pyrolytic graphite (HOPG) are 
obtained by cooling the sample to 1 IK and operating at large 
tip- sample distances — the first time such a measurement has 
been performed despite the long history of interest in the K- 
on-graphite system. The pair distribution function extracted 
from the spatial arrangement of K atoms deviates strongly 
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FIG. 1: Topography of the HOPG graphite surface prior to potassium 
deposition, [(a)] Large scale image showing the clean surface. |(b)| 
Atomic resolution image of the area highlighted in[(a)]The tunneling 
bias and current are 100 mV and 100 pA. White rectangles are 5 nm 
(a) and 1 nm (b) scale bars. 



from a random distribution, although no long-range order is 
observed. Quantitative comparison to molecular dynamics 
simulations suggests that K-on-graphite carry a dipole mo- 
ment of 10.5±1 Debye due to the charge transfer. The ex- 
pected hexagonal rotational order for such isotropic repulsive 
interaction is broken by thermally activated diffusion. Time 
lapse imaging and atom tracking techniques allow diffusion 
trajectories to be mapped out in real time. From the diffusion 
rate, an energy barrier of 32 meV is estimated for K hopping 
between lattice sites. 

Measurements were performed under ultra high vacuum 
conditions (« 10" 10 mbar) at 11 K. A sample of HOPG 
was cleaved in-situ to avoid any possible surface contamina- 
tion from exposition to ambient atmosphere. A beetle type 
STM provided scanned images, using a W tip prepared in-situ 
through argon sputtering and high temperature annealing. A 



typical surface obtained after cleaving can be seen in Fig. 1(a) 



large atomically flat terraces free of defects impurities or sur- 
face contaminant confirm the instrument stability and surface 
preparation protocol. Figure |T(b)| zooms in a smaller region 
with atomic resolution, showing the triangular lattice charac- 
teristic of Bernal stacking where only one out of two atoms is 
observed fT3l . 
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Potassium deposition was carried out with a 20 sec dose 
from an SAES getter heated with 6.5A of current and located 
50 cm from the sample. After K deposition, the tip-sample 
distance was greatly increased compared to the standard tun- 
nelling conditions used for Fig. [I] in order to image the potas- 
sium atoms without perturbing them with the tip (dragging or 
pushing them around) . This was accomplished by operating 
at a high bias voltage (typically -2 V) and very low current 
(typically 5 pA) — the equivalent of 400 Gft impedance com- 
pared to 1 Gft in Fig. [I] Under those extreme tunnelling con- 
ditions, single potassium atoms are clearly observed as bright 
spots on the graphite sample (Fig. [2|. The drawback of the 
large tip-sample distance is the loss of atomic resolution in 
the substrate image. Calculations and LEED experiments in- 
dicate that K adatoms sit on the empty sites in the middle of a 
graphite honeycomb lattice cell (T4lfT7l . but this could not be 
confirmed due to the lack of atomic resolution in the graphite 
while imaging K atoms. 



A linescan through Fig. 2(a) indicates an effective height 
for the K atoms of 2 A and a width of 7.8 A (full width 
at half maximum), reminiscent of single adatoms on metal- 
lic surfaces [18-20]. By comparison, calculations and LEED 
measurements of the crystallographic height (center of the K 
adsorbate to center of the C surface atoms) suggest an actual 
height of 2.7 A(l4j \TE\ [XT). The discrepancy between these 
values may be explained by the fact than an STM measure- 
ment convolves the crystallographic height with tunnelling 
density of states integrated over bias. This explanation would 
imply that the density of states is smaller on potassium atoms 
than on the graphite surface. 

Topographic images like those in Fig. |2(a)| pro vide informa- 
tion about K-K interactions. It has been shown that graphite's 
two-dimensional and semi-metallic character strongly affect 
the spatial arrangement of potassium atoms on the surface, 
especially when compared to a metallic substrate (2D . At 



first glance the distribution of atoms in Fig. |2(a)| may ap- 
pear random, but statistical analysis demonstrates repulsion 
between adatoms. This is shown in Fig. 1 3 (a) | using a two di- 



mensional autocorrelation of STM topographs, where the dark 
area around the center indicates repulsion between potassium 
adatoms. For larger separations a rapidly decaying oscilla- 
tory behavior is observed, similar to what one might expect 
for a liquid with no long range orientational order. This struc- 
ture is typical for particles subject to a long range repulsive 
interaction in two dimensions l22l . and consistent with the 
diffraction patterns observed in LEED experiments for potas- 
sium deposited on graphite 0. 

The strength of the repulsive interaction can be estimated 
by matching the measured pair distribution function to one 
generated by simulation based on a particular pair poten- 
tial. Figure |3(c)| shows the measured distribution function 
compared to the distribution for noninteracting atoms, and 
to distributions generated by a molecular dynamics simula- 
tion (LAMMPS) for K atoms at UK constrained to move on 
a 2D plane (23). The simulations were perfomed assuming 
a dipole-dipole interaction (Vkk oc r -3 ), reflecting charge 
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FIG. 2: | (a) | Topographic image of the graphite surface after the de- 
positon of potassium, showing isolated potassium atoms. The tun- 
neling bias and current are -2V and 5 pA. |(b)| The height profile is 
taken along the dotted line drawn in[(a)] |(c)]3D landscape (full height 
200pm) of a single potassium atom. 



transfer from K that is completely screened by an image 
charge in the graphite |24, 25 ]. The family of simulated dis- 



tributions in Fig. 3(c) represents a range of dipole moments, 
always fixed in orientation to point out of the plane (23). As 
expected, a larger dipole leads to a more organized system, 
with sharper peaks in the pair distribution function at multi- 
ples of the nearest neighbor separation. The best fit to the 
experimental distribution is obtained with a dipole of 10.5 ±1 
D. The data could also be fit with Coulomb interaction be- 
tween monopoles (Vkk oc r _1 ), but the dipole interaction is 
more appropriate for this adatom density (typical spacing 3 
nm) because the expected screening length in graphite/ene is 
only a few lattice spacings, less than 1 nm (24l[25). The exact 
form of the potential is nevertheless likely to be more compli- 
cated than a simple dipole, and some deviations between sim- 
ulation and experiment were observed. For example, the pair 
probability below the first peak in the distribution (around 2 
nm) was consistently higher in the experiment than the dipole 
simulation. Electron-mediated oscillatory interaction between 
adatoms, commonly observed on metallic surfaces (26U28) , 
are not expected here because graphite does not support the 
same kind of surface states. 

The value of the dipole moment for screened potassium that 
is extracted from these measurements can be compared to pre- 
vious experiments and calculations. Reference 10 measured 
the evolution of the work function as a function of the cov- 
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FIG. 3: Study of the spatial organization of the diluted phase of the 
potassium atoms on graphite. |(a)| Autocorrelation image C(x,y) av- 
eraged on about 30 images similar to figure |2(a)| The dark area 
around the central peak evidences the repulsive interaction between 
the potassium atoms. |(b)|radi al distribution function C(r), i.e average 
over the angle of C(x,y). |(c)| Pair distribution function, i.e., probabil- 
ity of having two atoms separated by a given distance. The empty cir- 
cles are the experimental data; dotted lines are molecular dynamics 
simulations for evenly spaced dipole moment {6.2D..16.3D}, with 
a solid line for 10.5D; dashed line is the expected distribution for a 
system without interaction, i.e. setting the dipole moment to in the 
simulations. The inset shows the data and the best fit on a larger scale 
to show the effect of the finite size of the image. 



erage and reported a value of 9.4±1.5D, in very good agree- 
ment with our estimate. Calculations predict a value of 8.3 
D at low coverage (as would be applicable here), decreas- 
ing down to 4.5 D at high coverage due to a smaller charge 
transferlfTTl [29l . EELS measurements suggest that 0.7e is 
transferred per adatom to graphite in the diluted phase, consis- 
tent with theoretical calculations predicting between 0.4 and 
0.7e(l6l [Hi [301. Using 2.7 A as the height of the K atom 
above the graphite surface, 3.4 A as the interlay er spacing, 
and a charge transfer of 0.7e, the dipole moment extracted 
from our measurements implies that the screening charge lies 
almost entirely in the uppermost plane lfT6l . 

The interaction between potassium adatoms defines their 
instantaneous spatial distribution statistics, but the dynamics 
are governed primarily by interactions with the underlying lat- 
tice. At UK, diffusion of K on the graphite surface was slow 
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FIG. 4: : Diffusion of potassium atoms on the graphite surface. 
[(a)] Initial topographic STM image and trajectory of several selected 
potassium adsorbates (for clarity) over 2 hours (solid lines). |(b)| Evo- 
lution of the standard deviation a(t) of the distribution as a function 
of time. The solid line is a fit with a(t) — V2Dt, used to extract D. 
[(c)] shows the maximum value of the distribut ion P (0,t) as a function 
of time with a fit with P(0, t) 
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enough that most of the atoms remained localized while the 
STM tip was scanning over them, but fast enough that their 
their movement could be tracked over a timescale of minutes. 



One notices, in Fig. |2(a)| that a small fraction of adatoms ap- 
pear to be cut in half: these are atoms that moved to a differ- 
ent site as they were being imaged. If a scan were repeated 
immediately after the previous frame had completed, it was 
observed that most of the adatoms had hopped to a new site 
during the course of the full scan frame. From this we can con- 
clude that diffusion is dominated by thermal activation, and is 
not a tip induced artifact. 

The diffusion process was characterized quantitatively by 
recording a time-lapse sequence of frames over two hours, 



with each frame acquired in a few minutes. Figure 4(a) shows 
the first (t = 0) of a sequence of scans with K atoms at posi- 
tions fi(t = 0), superimposed by trajectories, fi(t), of several 
of the atoms. Considering an ensemble of many such trajec- 
tories, one obtains a probability distribution of the displace- 
ment, P(Sr,t), with Sri(t) = \fi(t) - fi(t = 0)|, that can 
be compared to the expected distribution for two dimensional 
diffusion: 



P(5r,t) 
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a normalized gaussian in Sr with a time-dependent standard 
deviation a for diffusion constant D. The diffusion constant 
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was extracted from the data in two ways. First, the measured 
distribution at each time was fit to a gaussian, an d the r esulting 
standard deviation was fit to a(t) = \/2Dt (Fig. 4(b) ), giving 
L> = 2.6±0.8xl0 _4 nm 2 /s. Next, the experimentally mea- 
sured probability of the atom to remain at its initial position 
was fit to P(0, t) = er /( tJ^j) ' where erf is the error func- 



tion and ro the size of the bins used in the analysis (Fig. 4(c) ). 
This gives D = 2.3 ± 0.7 x 10 ~ 4 nm 2 /s. Combining both, we 
obtain D = 2.4 ± 0.5 x 10" 4 nm 2 /s. 

The diffusion coefficient can be related to the temperature, 
T, by a simple thermally activated hopping model PTI : 



D 



k B Tnl 2 ^jl 

e k B T 

2ha 



(2) 



where k# is the Boltzman constant, h is the Planck con- 
stant, n = 6 is the number of neighbouring lattice sites avail- 
able to a particle on graphite, a = 2 is the dimensional- 
ity of this system, E is the energy of the diffusion barrier 
and I = 2.46 A is the spacing between lattice sites, giv- 



ing D/T = 1.5 



W 10 e -E/ kl 



T 2 

nm s 



L K for graphite. 



The diffusion coefficient extracted above, D=2.4±0.5 x 10 -4 
nm 2 /s, at a measurement temperature T=11±1K, implies an 
energy barrier E = 32±3 meV (Fig. 4(d)) . By comparison, 
the interaction energy between dipole moments of 10. 5D at a 
distance of 3 nm is 2.5 meV. One can conclude from this that 
lowering the temperature would lead to an organized structure 
of potassium atoms only for coverages significantly higher 
than those reported here (32). This is the first experimental 
measurement of the hopping barrier for K on graphite, and is 
consistent with recent predictions of E « 50 meV (TSJ 02) 
based on density functional calculations. 

In conclusion we have reported the imaging of single K 
atoms adsorbed on HOPG by scanning tunnelling microscopy. 
We have used standard statistical methods to extract the po- 
tential landscape for a single adsorbate. It is composed of two 
main contributions: on the one hand an isolated potassium ad- 
sorbate will see the diffusion barrier of 32 meV on the HOPG 
surface; the presence of other adsorbates adds a dipole-dipole 
electrostatic interaction. By comparison with molecular dy- 
namics simulations, our data suggests a value 10. 5D for the 
screened dipole moment. This kind of local probe study could 
bring new insights into phase transitions occurring in two di- 
mensions such as the ones observed for potassium adsorbed 
on graphite. Furthermore, the possibility to study single alkali 
atoms on the surface of graphite should pave the way towards 
the study of the local modification of the electronic structure 
induced by alkali atoms, which can contribute to the under- 
standing of mechanisms driving the macroscopic phenomena 
observed in the graphite/alkali system. 
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